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Abstract 

We consider solutions to nonlinear hyperbolic systems of balance laws with stiff relaxation 
and formally derive a parabolic-type effective system describing the late-time asymptotics of 
these solutions. We show that many examples from continuous physics fall into our framework, 
including the Euler equations with (possibly nonlinear) friction. We then turn our attention 
to the discretization of these stiff problems and introduce a new finite volume scheme which 
preserves the late-time asymptotic regime. Importantly, our scheme requires only the classi- 
cal CFL (Courant, Friedrichs, Lewy) condition associated with the hyperbolic system under 
consideration, rather than the more restrictive, parabolic-type stability condition. 
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1 Introduction 



This short presentation is based on the joint work [3] in coUaboration with C. Berthon and R. 
Turpault. We are interested in hyperboUc models arising in continuum physics and, especiahy, 
describing complex multi-fluid flows involving several time-scales. The partial differential equations 
under consideration are nonlinear hyperbolic systems of balance laws with stiff relaxation sources. 
We investigate here the late-time behavior of entropy solutions. 

Precisely, our objective is to derive the relevant effective system — which turns out to be of 
parabolic type — and to investigate the role of a convex entropy associated with the given system 
of balance laws. As we show it, many examples from continuous physics fall into our framework. 
In addition, we investigate here the discretization of such problems, and we propose a new finite 
volume scheme which preserves the late-time asymptotic regime identified in the first part of this 
paper. 

An outline of this paper follows. In Section 2, we present a formal derivation of the effective 
equations associated with our problem. In Section 3, we demonstrate that many examples of 
continuous physics are covered by our theory. Finally, in Section 4, we are in a position to present 
the new discretization method and state several properties of interest. 

2 Late-time/stifF relaxation framework 

2.1 Hyperbolic systems of balance laws 

We consider systems of partial differential equations of the form 

edtU + da:F{U) ^ U ^ U{t,x) e n cR^ , (2.1) 

where t > 0, x G M denote the (time and space) independent variables. We make the following 
standard assumptions. The flux F : il is defined on a convex and open subset fi. The 

first-order part of (|2.ip is an hyperbolic system, that is, the matrix-valued held A{U) :— DuF{U) 
admits real eigenvalues and a basis of eigenvectors. 

We are interested in the singular limit problem £ — >■ in the limit of late-time and stiff 
relaxation. In fact, two distinct regimes for systems like (|2.1|) can be considered. In the hyperbolic- 
to-hyperbolic regime, one replace edtU by dtU , and one establishes that solutions to 

dtU + d,F{u) = -^^^, u^uit,x), 

e 

are driven by an effective system of equations (£ — )■ 0) of hyperbolic type. Such a study was 
pioneered by Chen, Levermore, and Liu [5]. On the other hand, in the hyperbolic-to-parabolic 
regime under consideration in the present work, we obtain effective equations of parabolic type. 
Earlier work by Marcati et al. discussed this regime too and established several important 
convergence theorems. 

Our objective here is to introduce a general framework to deal with such problems. We make 
the following assumptions. 

Assumption 1. There exists an n x iV matrix Q with (maximal) rank n < N such that 

QRiU) ==0, u en. 

Hence, QU e Q^, —: uj satisfies 

edt{QU) +d,{QF{U)) =0. 
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Assumption 2. There exists a map £ : u) C — )• f2 describing the equilibria u € w, with 

R{£{u))=0, u = Q£{u). 

It is then convenient to introduce the equihbrium submanifold M. := {J7 = f (u)}. 

Assumption 3. It is also assumed that 

QF{£{u)) =0, ueuj. 

To motivate this condition we observe that, formaly at least, the term dx{QF{£{u))^ must vanish 
identitically, so that QF(£{u)) must be a constant, conveniently normalized to be 0. 

Assumption 4. For all u G co, 

dim(^keT{B{£{u)))j = n, 
kcr(B(£(u))) nim {B{£{u))) = {0}. 
Hence, the N x N matrix B := DRu has "maximal" kernel on the equilibrium manifold. 

2.2 Chapman-Engskog-type expansion 

We proceed by using a Chapman-Engskog expansion in order to derive effective equations satisfied 
by the local equilibria u = u{t, x) G u. So, we write 

=£iu)+eUi+e''U2 + ..., u:=QW, 

which should satisfy e 5* f/^ + d^F{U^) = -R{U^)/e. It follows that 

QUi=QU2 = ... = 0. 

For the flux we find 

F{U') = F{£iu)) + e A{£{u)) Ui + 0{e''), 

and for the relaxation 

?^ = B{£{u)) C/i + ^-DlR{£{u)).{U^,U^) + eB{£{u)) + 0{e^). 
In turn, we deduce that 



edt{£{u))+dx{F{e{u)))+edx[A(£{u))U, 
= -B{£{u)) C/i - ^-DlR{£(u)).{lJuVi) - eB{£{u)) U2 + 0{e^). 

We begin by considering the zero-order terms and thus deduce that Ui G satisfies the linear 
system 

B{£{u)) f/i = -54F(f («))) e M^. 
We can solve this equation in Ui provided we recall that QUi = and observe the following. 

Lemma 2.1 (Technical lemma). Let C be an N x N matrix satisfying dimker C = n, and ker Cfl 
ImC = {0}, and let Q be an nx N matrix of rank n. Then, for all J G K^, the system 

CV = J, 
QV = 0, 

admits a unique solution V G M™ if and only if QJ = 0. 
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We can thus conclude with the following. 

Proposition 2.2 (First-order corrector problem). The first-order corrector Ui is characterized 
uniquely by 

B{S{u))Ui = -d4F{S{u))), 
QUi = 0. 

We now turn our attention to the first-order terms and we arrive at 

dt{£{u)) +d^(^A{£{u))Ui) = -^DlR{£{u)).{Ui,Ui) - B{£{u})U2. 

Multiplying by Q and using Q£(u) = uwe find 

dtu + d^(QA{£{u)) Ui) = -^QDfjR{£{u)).iUi, f/i) -QB{£{u)) U2. 

But, by diff'erentiating the identity QR{U) = 0, it follows that 

QDfjR.{Ui,Ui) = 0, QBU2 = 0. 
Theorem 2.3 (Late time/stifi^ relaxation effective equations). One has 

dtU = -d^(QA{£{u)) C/i) =: {M{u)d^u) 

for some nx n matrix M{u), where U\ is the unique solution to the first-order corrector problem 

B{£{u))Ui = -A{£{u))d^{£{u)), 
QUi = 0. 

2.3 The role of a mathematical entropy 

Next, we investigate the consequences of assuming the existence of a mathematical entropy $ : 
— )■ K, satisfying by definition: 

Assumption 5. There exists an entropy- flux \I> : — ^ R such that 

Du^A = Du'if inO. 

So, all smooth solutions satisfy 

eat($(C/^)) +5,(*(t/^)) = -Du-fiU') 

and, consequently, the matrix D^<& A is symmetric in fl. In addition, we assume that the map $ 
is convex, i.e. the N x N matrix Df^^ is positive definite on Ai. 

Assumption 6. The following compatibility property with the relaxation term holds: 

Du^R>0 infi, 

Du^U) = v{U)Q e K^, v{U) e W^. 

Returning to the effective equations 

dtU = dj), V := -QA{£{u)) Ui 

and multiplying it by the Hessian of the entropy, we conclude that the term Ui G is now 
characterized by 

£{u)U, = -{Dl<P){£{uj)d4F{£{u))), 
QUi = 0, 
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where C{u) = Dl^{£{u))BiS{u)). 

Then, using the notation £{u)~^ for the generaUzed inverse with constraint and setting 



we obtain 

Finally, one can check that 



with V := {Du<i>{S) f. 



S{u) := QA{S{u)), 
V = S£-'{Dl<!>){£)d4F{£)). 



Theorem 2.4 (Entropy structure of the effective system). When the system of balance laws is 
endowed with a mathematical entropy, the effective equations take the form 

dtU = d,(^L{u)d4D^<l>i£{u))f), 

where 



Liu) 



S{u)C{u)-^Siu) 



T 



= QA{£{u)), 

= {Dl^){£{u))B{£iu)), 



S{u) 
L{u) 

where, for all b satisfying Qb = 0, the unique solution to 

jC{u)V = b, QV = 

is denoted by C{u)^^b (generalized inverse). 

Alternatively, the above result can be reformulated in terms of the so-called entropy variable 
(^Du^{£{u))) . Furthermore, an important dissipation property can be deduced from our as- 
sumptions, as follows. Prom the entropy property and the equilibrium property R{£{u)) = 0, we 
find 

Du^R > inf2, 
(Du^R) \u=siu) =0 in w. 

Thus, the matrix Dfj 

u=s(u) is non-negative definite. It follows that 
Dl [Du^R) = Dl^B + [Dl^Bf when U = £{u), 

so that 

Dl^B\u^e(u)>^ inw. 

The equilibrium entropy ^{£{u)) has the property that its associated (entropy) flux u 
^{£{u)) is constant on the equilibrium manifold oj. Indeed, for the map '^{£), we have 

Du{^{£)) = Du^{£)Du£ 

= DuH£)Ai£)Du£. 

Observing that {Du<i>){£) = Du{^{£))Q, we obtain 

D„ (*(£(«))) = D^^£{u))QA{£{u))Du£{u) 
= Du{^{£{u)))d^QF{£{u)). 

Since QF{£) = 0, then DuQF{£) = and the proof is completed. 
Therefore, £>„(*(£ (u))) = for all u e w. 
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Recalling the asymptotic expansion 

W = S{u)+£Ui + 
where Ui is given by the first-order corrector problem, we write 

and deduce 

d,,^{U^')=ed^^Du-^{,S{u))Ui+0{e'). 
Similarly, for the relaxation source term, wo have 

Du^{U')R{U^) = e^Dl^{£{u))DuR{£{u))Ui + 0{e^). 

We thus obtain 

dt{^{£{u))) + d,(Du-^{£{u)) Ui) 
= -U^ {Dl^{£{u))B{£{u))) C/i. 
But, we have already established 

X {Dl^){£)B{£)X>Q, XeK^. 
Proposition 2.5 (Monotonicity of the entropy). The entropy is non-increasing in the sense that 

dtm£{u))) +d,iDu^{£{u))Ui)<0. 
In the notation given earlier, one thus have 

dtm£{u))) = d,(^{D4^£{u)))L{u)d,{Du{^{£{u))f). 

Remark 2.6. The above framework can be extended to handle certain nonlinear diffusion regime, 
corresponding to the scaling 

sdtU + d,F{U) = -^. 

The parameter q > 1 introduces an additional scale, and is indeed necessary for certain problems 
where the relaxation is nonlinear. The relaxation term is supposed to be such that 

R{£{u)+eU) =e'^R{£{u)+M{e)u), U € CI, ueoj, 

for some matrix M(e). In that regime, the effective equations turn out to be o/ nonlinear parabolic 
type. 

3 Examples from continuum physics 
3.1 Euler equations with friction term 

The simplest example of interest is provided by the Euler equations of compressible fluids with 
friction: 

edtp + dx{pv) = 

edtipv) +dx{pv^ +p{p)) = ^ 

in which the density p > and the velocity component v represent the main unknowns while 
the pressure p : R+ R+ is prescribed and satisfy the hyperbolicity condition p'{p) > for all 
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p > 0. Then, the first-order homogeneous system is strictly hyperbolic and p.ip fits into our 
late-time/stiff relaxation framework provided we set 

and Q = (1 0). In this case, the local equilibria u = p are scalar- valued, with 

and we do have QF{£{u)) = 0. 

The diffusive regime for the Euler equations with friction is analyzed as follows. First, according 
to the general theory, equilibrium solutions satisfy 

atp--a.(QA(£(«))c/i), 

where 

DuF{£{u)) = 

Here, Ui is the unique solution to 



1 

p'{p) 



B{E{u))Ui = -d^{F{E{u))), 
QUi^O 

with 

' \ , /^,^, ,A f 



B{£iu))=(^ ? )' ^-(^(^(^))) = ( 



dx{pip)) 

This leads us to the effective diffusion equation for the Euler equations with friction: 

dtp^dl{p{p)), (3.2) 

which is a nonlinear parabolic equation, at least away from vacuum, since p'{p) > by assumption. 
Interestingly, at vacuum, this equation may be degenerate since p'{p) typically vanishes at p = 0. 
For instance, in the case of polytropic gases p{p) = np'^ with k > and 7 G (1, 7) we obtain 

dtp = dx{p^^^dp). (3.3) 

In addition, by defining the nternal energy e{p) > by 

^ pip) 

e [P) = 

p 

we easily check that all smooth solutions to (13. ip satisfy 

edt[p^+peip)) +d,(^p'^ + ipeip)+p{p))v) = (3.4) 

so that the function 

HU) ^ p— + pe{p) 

is convex entropy compatible with the relaxation. All the conditions of the general framework are 
therefore satisfied by the Euler equations with friction. 
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3.2 Ml model of radiative transfer 

Our next model of interest arises in the theory of radiative transfer, i.e. 

edte + d^f = , 

£ 

/ 



edtf + d,{x{f/e)e)=-^, (3.5) 



e 

sdtT = , 

where the radiative energy e > and the radiative flux / are the main unknowns, restricted by 
the condition |//e| < 1, while r > denotes the temperature. The function x : [— 1, 1] — )• M+ is 
called the Eddington factor and, typically. 




5 + 2^4 - 3^2 • 

This system fits into our late-time/stiff relaxation framework if we introduce 

F{U) = Xipe j , R{U) J 
Now, the equilibria are given by u = r + r^, with 

£{u) =^Q^, Q := (1 1). 

We do have QF{£{u)) = and the assumptions in Sections 2 are satisfied. 
We determine the diffusive regime for the Ml model from the expression 






1 




( 


1 




x(o) 


x'(o) 


:h 







• 










^ 








where f/i is the solution to the linear problem 

C/i = 
(1 l)Ui = 0. 

Therefore, we have 






and the effective diffusion equation for the Ml system reads 



a,(T + T*)=a, (^T^a.r). (3.6) 



Again, an entropy can be associated to this model. 
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3.3 Coupled Euler/Ml model 

By combining the previous two examples we arrive at a more involved model: 

edtp + dx (pv) = 0, 
edtpv + d^{pv^ + p{p)) = --^P'" + 



(3.7) 



in which the same notation as before is used and k and a are positive constants. In the applications, 
a typical choice for the pressure is 

p{p) = Cpp", Cp < 1, 77 > 1. 

To fit this model within the late-time/stiff relaxation framework, we need to set 



U 



pv 
e 



( 



pv 



F{U) 



\ 



pv^ + p{p) 
f 

V X(|)e J 



R{U) 



( \ 

Kpv — af 



V '^f J 



The local equilibria are given by 







e 



u = QU ^ 






















I 


2 



and once again, one hSisQF{£{u)) = 0. 
We can then compute 



/ 1 0\ 

p'{p) 

1 

V i 0/ 







]:{-dMp)-\d..e) 



DuF{E{u)) = 

and we arrive at the effective diffusion system for the coupled Euler/Ml model: 



V 





3cr 



d^p--dlp{p)-^dle = Q, 
dte - ^dle = 0. 

OCT 



(3.8) 



The second equation is a standard heat equation, and its solution serves as a source-term in the 
first equation. 

3.4 Shallow water with nonlinear friction 

Our final examples requires the more general theory of nonlinear relaxation mentioned in He- 
mark[n?ni smd reads 

edth + dx{hv) =0, 

/ \ K'^(h) (3-9) 
edt{hv)+djhv^ +p{h)) = ^^ghv\hv\, 
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where h denotes the fluid height and v the fluid velocity v. The pressure is taken to be p{h) — 
gh^ /2. and 5 > is called the gravity constant. The friction coefficient k : M+ M+ is a positive 
function, and a standard choice is = ^ with kq > 0. 

The nonlinear version of the late-time/stiff relaxation framework applies if we set 

hv)' ^^^^ = ( hv^ + p{h) ) ' ^^^^ ^ ( K^h)ghv\hv\ 
The scalar equilibria u = h are associated with 

£iu) = ( ) , g = (1 0). 

Here, the relaxation is nonlinear and satisfies 

R{£{u) + eU) = e^R{£{U) + M{e)U) , 

with 

in turn , we may derive a nonlinear effective equation for the Euler equations with nonlinear 
friction, that is, 

d,h ^dj^ -4^] , (3.10) 

which is a nonlinear parabolic equation. 

In addition, by introducing the internal energy e{h) :~ gh/2, we observe that all smooth 
solutions to p.9p satisfy the entropy inequality 

edt (/.y + ffy) + 9. [h^ + gh^^ V = -^ghv^\hv\. (3.11) 

The entropy 

satisfies the compatibility properties relevant to the nonlinear late-time/stiff relaxation theory, 
with in particular 

i?(£:(.) + M(0)C7,)^( 
where Ui = (0 ^) . We obtain here R{£{u) + M(0)f7i) = c{u)Ui with 



c(u) — gK{h)\Jh\dxh\ > 0. 



4 Asymptotic-preserving finite volume schemes 
4.1 General strategy 

We are going now to design finite volume schemes that are consistent with the asymptotic regime 
e — > determined in the previous section and, indeed, to recover an effective diffusion equation 
that is independent of the mesh-size. The discretization of hyperbolic-to-hyperbolic regimes were 
investigated first by Ji and Xin [12^. Here we propose a framework to cover hyperbolic-to-parabolic 
regimes. For earlier work on this latter regime see [U [21 El [7] . 

Step 1. Our construction is based on a standard finite volume scheme for the homogeneous system 

dtU + d.,F{U) = 0, 

and we will begin by describing such a scheme. 
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Step 2. We then will modify the above scheme and include a matrix-valued free parameter, allowing 
us to approximate the non-homogeneous system 

dtU + d^F{U) = --fRiU), 

for a fix coefficient 7 > 0. 

Step 3. Finally, we will perform an asymptotic analysis after replacing the discretization parameter 
At by eAt, and 7 by 1/e. Our goal then will be to determine the free parameters so that to 
ensure the asymptotic-preserving property. 

Let us briefly define the so-called HLL discretization of the homogeneous system, as proposed 
by Harten, Lax, and van Leer [11]. For simplicity we present here the solver based on a single 
constant intermediate state. The mesh is assumed to be uniform mesh made of cells of length Ax: 

r 1 

R-i/2,a;j+i/2j, = Xi + — 

for all? = ... , —1, 0,1,.. .. The time discretization is based on a parameter At, restricted by the 
famous CFL condition (Courant, Friedrichs, Lewy, cf. 4 ) with 

Starting from som initial data (lying in the convex set fl): 

we construct piecewise constant approximations at each time i™: 

= Ur, xe [a;.- 1/2, 2^^+1/2), I e Z. 

Following Harten, Lax and van Leer [llj , at each cell interface we use the approximate Ricmann 
solver: 

Ul, t < "b, 



t 



Ur, |>5, 

where 6 > is sufficiently large. The "numerical cone" (and numerical diffusion) is thus determined 
by the parameter 6 > and here, for simplicity in the presentation, we have assumed a single 
constant b but, more generally, one can introduce two speeds b~^^^2 < ^i+1/2 each interface. 
We introduce the intermediate state 

U* - 1{Ul + Ur) - Y^{FiUR) - F{Ul)), 

and assume the CFL condition 

^A< ^ 1 
- 2' 

so that the underlying approximate Riemann solutions are non-interacting. Our global approxi- 
mate solutions 

UAxix^t"" + t), ie[0,At), xeM. 

are then obtained as follows. 

The approximations at the next time level i™"*"^ are determined from 

ur^' = ^ r^''' UAAx,t"^ + ^t)dx. 



-1/2 
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Then, recalling i7*+i/2 = M^r + ^i+i)- 5E(^(^i+i)--P'(f^r)): and integrating out the expression 
given by the Riemann solutions, we arrive at the scheme for the homogeneous system 

TTm+l jjm At / j-,HLL jpHLL \ 

where 

K"i/2 = I {nun + nuiii)) - \(ur+, - un- 

(More generally one could take into account two speeds bT'_^^^^ < 

We observe that the above scheme enjoys invariant domains. The intermediate states U*^-^^2 
can be written in the form of a convex combination 

provided b is large enough. An alternative decomposition is given by 

where A is an "average" of DjjF. By induction, it follows that in fl for all m, i. 

4.2 Discretization of the relaxation term 

We start from the following modified Riemann solver: 

Ul, f<-6, 
U*^, -6 < I < 0, 

< I < 6, 
Ur, I > 6, 
with the following states at the interface: 

U*^ = qU* + {I-a) {Ul - R{Ul)) , 
U*"^ = aU* + {I-a) [Ur - R{Ur)) . 
We have here introduced some N x A''-matrix and N-vectov defined by 



Uni^;UL,UR) = { 



+ \ R{U) = {I + a)-'R{U). 



The term ct is a parameter matrix to be chosen so that (all inverse matrices are well-defined and) 

the correct asymptotic regime is recovered at the discrte level. 
At each interface we use the modified Riemann solver 

and we superimpose non-interacting Riemann solutions 

U^^{x,t"' + t), te[0,At), xeR. 

The new approximate solution at the next time f^+i is 



Ur^'= UZ{x,t"' + At)dx. 
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By integrating out the Riemann solutions, we arrive at the discretized balance law 

- ^{i - ^i+i/2)Ri+i/2iun- 

Observing that the discretized source can rewritten as 

^(/ - a,^,/,)Ri+i/2{un = ^«i+i/2(«r+i/2 - m+i/2{un 

= l<^^+i/2R{ur) 

and, similarly, 

we conclude that the proposed finite volume schemes for late-time/stiS-relaxation problems takes 
the form 

= ^{a,+i/2-^i-i/2)Fiun - 5fe+i/2 +«.-i/2)i?(t/r). 

Theorem 4.1 (Properties of the finite volume scheme). Provided 

Q-i+\l2-3.i-XI2 = 0{i\x) 

and the matrix-valued map a is sufficiently smooth, then the modified finite volume scheme is 
consistent with the hyperbolic system, with relaxation. The following invariant domain property 
holds: provided all the states at the interfaces 

utii/2 = »i+i/2ut+i/2 + {1- ^i+i/2){ur - mn), 
u:^i/2 = «i+i/2f/r+i/2 + {1- m+i/2){ur+, - mr+i)) 

belong to Vt, then all the states f/|" belong to O. 

4.3 Discrete late-time asymptotic regime 

As explained earlier, we now replace At by At/e and 7 by 1/e, and consider 

= ■^i^i+i/2-^i-i/2)Fiun - ^(«i+i/2 +«.-i/2)^(t/r), 



where 



Ax ^ ~^ 

QLi+1/2 = (-'" + + ^i+1/2) 



2er 

Plugging in an expansion near an equilibrium state 

Ur =S{u^)+e{U{)T + 0{e''), 
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we find 

The first-order terms for the discrete scheme leads us to 



2b 



(^{1 + 0:2+1/2) ^-Pi+i/2l£(") - (^ + ^1-1/2) ^-Pi-i/2l£(«) 
+ ^ ((^ + -{1 + ^i-1/2)-') F{£{<)) 

- ^ {{I + a.i+,/2)-' + (/ + ^i-i/2)-^) B{£{uT))mT. 
At this juncture, we assume the existence of an n x n matrix such that 

We then multiply the equation above by the n x N matrix Q and obtain 



At 

with 



ur ■'-ur) = (Mi+y2QF,1if2\siu) - Mi.y^QFlLkfMu)) 



The discrete asymptotic system is thiis 

- O = -l^(M,+y2{nT+i - ur)+M,.y2{uT-i ' <))• (4-3) 
Recall that for some n x n matrix M (u) the effective equation reads 

dtu = dx {M{u)dxu) . 

Theorem 4.2 (Discrete late-time asymptotic-preserving property). Assume the following condi- 
tions on the matrix-valued coefficients: 

• The matrices ^ 



2 



are invertible for e G [0, 1]. 

• There exists annx n matrix A1,+i/2 ensuring the commutation condition 

• The discrete form of M.{u) at each cell interface a;i+i/2 satisfies 

Mi+x/2 = M{u) + 0{Ax). 
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Then, the effective system associated with the discrete scheme coincides with the one of the late- 
time/stiff relaxation framework. 

We refer the reader to [3] for numerical experiments with this scheme, which turns out to 
efficiently compute the late-time behavior of solutions. It is observed therein that asymptotic 
solutions may have large gradients but are actually regular. We also emphasize that we rely here 
on the CFL stability condition based on the homogeneous hyperbolic system, i.e., a restriction on 
At/ Ax only is imposed. In typical tests, about 10000 time-steps were used to reach the late-time 
behavior and, for simplicity, the initial data were taken in the image of Q. A reference solution, 
needed for a comparison, was obtained by solving the parabolic equation, under a (stronger) 
restriction on At/{Ax)'^. 

The proposed theoretical framework for late-time/stiff relaxation problems thus led us to the 
development of a good strategy to design asymptotic-preserving schemes involving matrix-valued 
parameter. The convergence analysis (e — 0) and the numerical analysis (Aa; — 0) for the 
problems under consideration are important and challenging open problems. It will also very 
interesting to apply our technique to, for instance, plasma mixtures in a multi-dimensional setting. 

Acknowledgments 

The author was partially supported by the Agence Nationale de la Recherche (ANR) through the 
grant 06-2-134423, and by the Centre National de la Recherche Scientifique (CNRS). 

References 

[1] C. Berthon, P. Charrier, and B. Dubroca, An HLLC scheme to solve the Ml Model of radiative 
transfer in two space dimensions, J. Sci. Comput. 31 (2007), 347-389. 

[2] C. Berthon and R. Turpault, Asymptotic preserving HLL schemes, Numer. Methods Partial Differ- 
ential Equations, DOL 10.1002/num.20586. 

[3] C. Berthon, P.G. LeFloch, and R. Turpault, Late-time relaxation limits of nonlinear hyperbolic 
systems. A general framework. Preprint arXiv:1011.3366 

[4] R. Courant, K. Friedrichs and H. Lewy, Uber die partiellen Differenzengleichungen der mathematis- 
chen Physik, Math. Annalen 100 (1928), 32-74. 

[5] F. Bouchut, H. Ounaissa, and B. Perthame, Upwinding of the source term at interfaces for Euler 
equations with high friction, J. Comput. Math. Appl. 53 (2007), 361-375. 

[6] C. Buet and S. Cordier, An asymptotic preserving scheme for hydrodynamics radiative transfer 
models: numerics for radiative transfer, Numer. Math. 108, 199-221 (2007). 

[7] C. Buet and B. Despres, Asymptotic preserving and positive schemes for radiation hydrodynamics, 
J. Comput. Phys. 215 (2006), 717-740. 

[8] G.Q. Chen, CD. Levermore, T.P. Liu, Hyperbolic conservation laws with stiff relaxation terms and 
entropy, Comm. Pure Appl. Math. 47, 787-830 (1995). 

[9] D. Donatelli and P. Marcati, Convergence of singular limits for multi-D semilinear hyperbolic systems 
to parabolic systems. Trans, of AMS 356, 2093-2121 (2004). 

[10] P. Marcati and B. Rubino, Hyperbolic to parabolic relaxation theory for quasilinear first order sys- 
tems, J. Differential Equations 162, 359-399 (2000). 

[11] A. Harten, P.D. Lax, and B. van Leer, On upstream differencing and Godunov-type schemes for 
hyperbolic conservation laws, SIAM Review, 25 (1983), 35-61. 

[12] S. Jin and Z. Xin, The relaxation scheme for systems of conservation laws in arbitrary space dimension. 
Comm. Pure Appl. Math. 45 (1995), 235-276. 

[13] P. Marcati, Approximate solutions to conservation laws via convective parabolic equations. Comm. 
Partial Differential Equ. 13 (1988), 321-344. 

[14] P. Marcati and A. Milani, The one-dimensional Darcy's law as the limit of a compressible Euler flow, 
J. Differential Equ. 84 (1990), 129-146. 



15 



